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ABSTRACT 

The hulk properties of tropical cloud clusters, such as the vertical mass flux, thd excess temperature, and> 
moisture and the. liquid water content, of the clouds, are determined from a combination of t the observed 
large-scale heat and moisture budgets over an area covering the cloud cluster, and 2) a model of a cumulus 
' ensemble which exchanges mass), heat, water vapor and' liquid water with the environment through entrain- 
ment and detrainment. The method also provides aft understanding of how the environmental air is heated 
"and moistened by the cumulus convection. An estimate of the average cloud cluster properties and the 
reheat and moisture balance of the environment, obtained from 1956 Marshall Islands data, is presented. 
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1. Introduction 

Advances in our understanding of the dynamics of 
tropical disturbances call for the clarification of a long- 
standing problem, i.e., the relationship between the ^ 
.large-scale motion and organized cumulus convecliof^ 
The important role played by deep cumulus clouds in 
the heat balance of the tropical atmosphere was first 
pointed out by Riehl and Malleus (1958). They showed 
that deep cumulus clouds can carry the released latent 
heat of condensation to the upper troposphere where 
. ; the local vertical gradient of moist static energy 
\Qt—c v T'\-gz-\-Lq) is positive. The role of penetrative 
■ - cumulus towers in the formation and growth of tropical 
*£y clones was discussed by Riehl and Malkus (1961) 
and by Yanai (1961a, b). They proposed that cumulus 
convection, as h ' heating mechanism, must be ad- 
equately parameterized in the framework of the large- 
-scale motion., The first theoretical > formulation of the 
cumulus, macro-motion coupling through boundary 
layer convergence was made by Charney and’Eliassen 
(1964) and Ooyama (1964, 1969a), in their attempts to 
explain the growth of tropical cyclones. . 

tin the case of weaker tropical disturbances, such as 
easterly waves (Riehl, 1945, 1948, 1954; Palmer, 1952), 
the connection with cumulus convection remained 
somewhat vague, although the control of cloud organiza- 
tion and precipitation by the large-scale motion had 
been suggested in a number of studies (e.g., Riehl, 
1954; Malkus, 1958; Malkus and Riehl, 1964). 

Based on an inspection of cloud, pictures taken by 
satellites during 1967, the JOC Study Group on Tropical 
Disturbances (1970) documented unquestionable evi- 
dence; of the organization of maritime tropical clouds 
into “cloud clusters.” The high degree of organization 
of '.colli vective clouds into clusters Suggested their 
likely association with large-scale Wave disturbances. 
Evidence of a systematic westward-propagation of the 
cloud clusters was presented bv Chang (1970). - - 
In parallel with the recognition of the tropical cloud 
clusters, rapid progress has been made in observational 
and theoretical studies of the tropical wave distur- 
bances [see review articles of Wallace (1971) and Yanai 
(197 la) J. We now have several working hypotheses 
concerning the origin of the waves and their interaction 
with organized cumulus convection. The so-called 
“wave, cloud-cluster interaction” has become a leading 
problem in tropical meteorology. 

Recent studies of the energetics of the tropical wave 
disturbances favor the view that the wave disturbances 
derive their kinetic energy primarily from the eddy- 
available potential energy which is generated by the 
release of the latent heat of condensation (Manabe 
and Sraagorinsky, 1967; Manabe el aL, 1970; Nitta, 
1970, 1972a; Wallace, 1971). The theory of how tropical i 
waves respond to the heating of the condensation has 
been advanced bv Qovama (1969b), Yamasaki (1969), 
Hayashi (1970), Holton (1971), T. Murakami (1972), 
and others. Holton et aL (1971) and Yamasaki (1971) 
further showed that the waves can be enhanced at 
•certain critical latitudes where the Coriolis parameter 
-eq uals the frequency of the waves. This is because of 
^ theYtngular response of the boundary layer which con- 
trols the large-scale moisture convergence. 



To further advance the theory of tropical waves, and 
to improve the numerical simulation of the tropicail 
general circulation, a much better knowledge of the 
properties of the tropical cloud clusters is needed. 
Recent!}', Williams (1970) and Reed and Recker (1971) 
studied the average macro -structure of the cloud 
clusters by making composites of many clusters. Reed 
'and Recker (1971) and Nitta (1972a) obtained vertical 
profiles of the heat source associated with the average 
cloud cluster. 

This study uses a somewhat different approach to 
obtain the average structure of a cloud cluster and its 
interaction with the environment. The study was made 
in parallel with Prof. Akio Arakawa J s effort to improve 
the parameterization of cumulus convection in the 
UCLA General Circulation Model (Arakawa, 1969, 
1971, 1972; Arakawa and Schubert, 1973). It is evident 
that a large amount of latent heat is liberated in cum- 
ulus clouds and the released heat is transported up- 
ward, but hma this' heal is utilized in warming the large - 
scale environment is not so obvious. We must have a 
better understanding of this physical process to im- 
prove cumulus parameterization. 1 

1 Earlier attempts of parameterization of cumulus convectiop^ 
were reviewed by Yanai (1971a), Bates (197.2) and Ogura (1970). 

Wc shall demonstrate that the observed large-scale 
heat and moisture budgets, together with a simple 
model of a cumulus ensemble, make it possible to 
obtain a large amount ;pf .information such as the, mass 
flux, the excess temperature, and the moisture and the 
liquid water content of the cumulus clouds. We shall ; 
then show how the observed heat and moisture. budgets' 
of the large-scale environment can be explained in 
terms of induced downward motion between the clouds' 
and the detrainment of water vapor and liquid water 
from the clouds. ' ‘ 

2. Method < ’’ ■ 


a. Preliminary considerations • 

For convenience, we define the dry static, energy as 
s^c'jX-t-gz, 'i ; (0 

and the moist sialic energy as ; 

li=CpT-\-gz-\-Lq, ■ ; (2) 


where 5 is the sum of enthalpy 1 and potential eneVgy'ahd 
k the sum of s and latent energy. For dry adkifcatic 


processes, 


j t ■ , 



while for both dry and moist adiabatic processes, 



dt 


We consider an ensemble of cumulus clouds which is 
embedded in a tropical large-scale motion system. 
Then we imagine a horizontal area that is large enough 
to contain the ensemble of clouds, but small enough to 
be regarded as a fraction of the large-scale system.- 
We now write the equations of mass continuity, heat 
energy and moisture continuity, averaged over .the 
hypothesized area: 

f- 'i'U 

V-V+— =0, . ... (5) 

dp ■/. 
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ds dso> 

— -f-V‘jrV4 — Qk^L{c— e), . (6) 

dt dp 

dq dgu 

bV-?V-f— = e-c. .(7) 

dt dp 

In the above, most of the notations are conventional; 
g> is the average vertical p velocity, Qn the heating 
rate due to radiation, c the rate of condensation per 
.-.unit mass of air, and e the rate of re- evaporation of 
cloud droplets. The horizontal averages are denoted by 

O- 

Eqs. (6) and (7) may be rearranged to give 

ds dSoi d 

^=__ +v *iV+— =Q R +L{c-e) (8) 

dt dp dp 

( dq — dqw\ d 

— t-V-tfVd— W(o'-e)+X~?V, (9) 

dl dp / dp 

where deviations from the horizontal averages are de- 
noted bv primes. It will be assumed that small-scale 
,<eddies in the horizontal components of wind have no 
-significant correlations with 5 ' and q r . But the presence 
■ of cumulus convection, and turbulent motion in the 
r.sub-cloud layer, may cause significant vertical eddy 
transports of heat and moisture. Eq. (8) shows that the 
apparent heating Q\ of the large-scale motion system 
consists of the heating due to radiation, the release of 
latent heat by net condensation, and vertical converg- 
ence of the vertical eddy transport of sensible heat. 

Eq. (9) is the equation of moisture continuity expressed 
in units of heating rate, with () 2 being a measure of the 
apparent moisture sink which is due to the net condensa- 
tion and vertical divergence of the vertical eddy trans- 
port of moisture. 

From (8) and (9), we obtain 

d d 

Qi~'Q 2 ~Qk— — ' (s'-j-Lq )&/ — — —ii<J , ( 10 ) 

dp dp 

where h! 0 / is a measure of the vertical eddy transport 
of total heat and may be used to measure the activity 
.of cumulus convection. For later convenience we define 
the vertical eddy flux of total heat by ^ 

F = [ P (Qt-Qz-Q^dp, (11) 

£ £ d vtp) - — — r 


where pr is the pressure at the cloud top. 

Integrating (8) from pr to po (the pressure at the 
j( .sea surface), — - 




(< Qi-Qx)dp 




(c e)dp (5 0} 



gJpfj £ 

=LP 0 +Sc, 


where Po and 5o are, respectively, the amount of pre- 
cipitation and the supply of sensible heal from the ocean 

<mrf'v/-A — — 




5 ~ (a fj.r f 


where Eo is the rate of evaporation from the sea. 

From (11), (12) and (13), we find _____ 

F 0 ^ F{p 0 ) = S*+LEf^^ (14) 


Eqs. (12)— (14) can be used to check the accuracy of 
the estimates of Q 1} @ 2 and Qr against surface observa- 
tions of Fa, So and E 0 . 



b. Cloud mass flux 

The large-scale vertical motion u> is the average of 
the mean ascent in active cumulus clouds and the slow 
descent or ascent of the environmental air. For sim- 
plicity we define the average mass flux across a unit 
horizontal area by 2 

s Dimensionally, it may be more logical to define the mass flux 
by M = —w/g ~pw. 

M = — w, (15) 

.and divide M into two parts, 

M=M c +tf, (16) 

M c being the upward mass dux in the active cumulus 
clouds, and M the residual mass flux in the environment. 
We assume that the active clouds occupy a fraction <r 
of the horizontal area. Then 


M c «= ~<XCi) c 

M=-( l-(r)« 




where u c is a typical vertical p velocity in the clouds 
and Ho a typical vertical velocity in the environment ; d) 
may be quite different from w. For example, if af 
ascent required by the large-scale conwu^«Tt(^takes- 
place in the cumulus clouds. <7 ®rimd §^0. If all 

the ascent is compensated by the descent between 
clouds, co=0 and (1 -a) d) = — crco c . 

We next assume that the dry static energy s and mix- 
ing ratio q have characteristic values inside the active 
clouds, which are denoted by and q c , respectively. 
The characteristic values of s and q in the environment 
are denoted by s and q. Therefore, 




S\- 

C*J 


S = aJ 5 -{-(l— <r)x 
q = oq t -\-{\ — <r)q 


( 18 ) 


It is important to note that 

I = a(s c — s) -J-Jss s 

tf=o-(} c — q)+q~q 

because 

(s c —s)<s and (q c ~ q)<q. 

Let us assume that above the cloud base the vertical 
eddy transports of s and q are solely due to the presence 
of active cumulus convection. Then, 



.s'c/ — so.' —sd: ~<r(l — tr) (s c — j) (w t — w) 
q'o)' ~qo>—qd)=cr(\ —a) (q c — q) (w c — w) 


( 20 ) 
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.Because of the foregoing reasons, we readily obtain 


—uuds c ~s) = M c {s c —s), (21a) 

o c (q c ~~q)-M c (q c —g). (21b) 

Eq. (21a) shows that the vertical eddy heat flux 
is given by the product of the cloud mass flux and 
■the temperature excess of the clouds, since s c — s 
~Cp(l\—T)\ Eq. (21b) has a similar interpretation. 
Such formulations were discussed by Riehl and Malkus 
(1961) and Yanai (1964). Combining (21a) with (2 lb), 
we have also 

-hW = s F=M c (h-h). (22) 

-Note that F is an observable quantity , as given by (U). 

The above derivations are correct only when all 
clouds at a given height share the same values of « c , 
s c and q c . In a more realistic situation, we may have a 
variety of clouds with different characteristics . We then 
generalize (21) and (22) (assuming ^pc 7 »<<Cl) by 

= ] C mfsci—s) 

___ i 

—q'u'=Y, ini(q ci —q ) 

-AV ntfKi-Jt), (24) 

i 

where i refers to the ith cloud. We note that 

Ic=E <»;. (25) 

1 

Thus, we now have 

d 

< Qi—Qr = ^12 {Ci—ei)+— 52 Mi(s c i~s), (26a) 

» dp i 

’ 3 

. Qt — L52 (d-ed-L—Jf mi(q c{ — q), (26b) 

» dp * 

where c,- and e 2 - are the rate of condensation taking place 
in the f th cloud^ind the rate of re-evaporation of that 
cloud, respectively. 




' X 

/>er €&yi-i}£‘ 


c, A model of the cumulus ensemble 

Consider an ensemble of cumulus clouds as shown in 
Fig. 1. These clouds are embedded in a large-scale 



Fig. 1. Schematic view of an ensemble of cumulus clouds. 

motion system. We shall assume that the activity of the 
cumulus clouds is statistically in an equilibrium with Ike 
imposed large-scale motion . Obviously, individual cumu- 
lus clouds will undergo their own life cycles. However, 
we hypothesize that, over a time scale describing the 
behavior of the large-scale motion system, the cloud 
ensemble maintains the heat and moisture balance with 
the environment. 

To simplify the treatment, we classify the cumulus 
clouds according to the heights of the cloud tops. 3 

3 This is equivalent to classifying the clouds according to their 
radii, because of the well-known entrainment-radius relationship. 

We assume that all the clouds have a common cloud 
base height, and that the values of s r and q r . (and thus 
k<), at a given height, are the same for every cloud 
which has the same top height. We may then use the 
suffix i as the index of a classified cloud type, instead of 
as the index of an individual cloud. The idealization thus 
introduced is shown in Fig. 2. 



Fig. 2. Idealized view of cumulus cloud types, classified according 
to their top heights (see text for notation). 


The equations of balance for mass, heat, water vapor 
and liquid water between cloud type i and its environ- 
ment arc written 



(mass) 


(heat) 


(water vapor) 


(liquid water) 


(27) 

(28) 

(29) 

(30) 


dp 
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where "SJ and are the rates of mass entrainment and 
detrainmen l per unit pressure interval. Because of the 
definition of the cloud mass flux, 8; and 5, have dimen- 
sion of the inverse of time ; c, is the rate of condensation, 
per unit pressure interval, occurring in cloud type i , U 
the mixing ratio of liquid water, and n the rate of gener- 
ation of raindrops which is assumed to fall out immedi- 
ately. In addition, we assume that the rate of re-evapor- 
ation ei is equal to that of the delrainmenL of liquid 
water , 4 

• 4 Thus the evaporation of falling raindrops is ignored. 

€; = dili . ( 31 ) 

The above set of. equations was presented by Arakawa 
(1971) and‘ Yanai (1971b). Ooyama (1971) also pre- 
sented similar equations. 

It is convenient to define the saturation moist static 
.energy of the environment as 

1i*=c p T+gz+I.q* (p,f), (32) 

where </* is the saturation mixing ratio at temperature 
T. Then it follows that 


(h ci -k*), (33) 

1+7 


L(q ni ~g*) = (hci-h*)> 

1+7 


(34) 


where 7 = (L/c r )(dq*/dT) p {hTter Ara.kawa(l969). 

We now have the unknowns m i} 5,-, 8,- ; s r .i, fhi> 
Ci, fi, When i M * x types of cumulus clouds are 

considered, we have (7X?mnx+2) equations £(27j (31), 
(33), (34) plus (26a, b)]. Obviously these are not 
sufficient to determine lOXimas: unknowns. 

In order to make the problem tractable, we introduce 
further simplifications, hirst we note that _ 

M c —Yl Si 




(35) 

(36) 


Then we define a weighted average of /z c ; by 

mikd/YL mjici/ M c , G^) 

i i i 

. and define 5 C , q e and i similarly. 

We now introduce the following assumptions: 


(i) Each type of cioud has a thin layer oj detrainmenl , 


where 

hd = h*. 


(38) 


This means that the clouds detrain at the level where 
they lose buoyancy. When (38) is inserted into (33), 
the criterion becomes 5 ,.,= 5, i.e., TV* = 7. It is also evi- 
dent from (34) that </,,- = y*. This assumption is also 
used by Arakawa (1971) and Arakawa and Schubert 
(1973). Near the tropopause we need a special con- 
. sideration because ibe cloud mass flux entering fiom 
below detrain in this layer regardless of the non- 
+u’oyancv condition (38). In this layer the entrai nmen t 
is assumed to be zero, so that 6~ dM c / d^TThls assump- 
tion allows the detrainment of heat as well as mois ture 
to take nl ace -m-th-is- lat^T 



(ii) The liquid water content of the detrained air, 
is the same as i.e., 

h=l (39) 

This gross assumption is needed to close the set of 
equations. 

(iii) The conversion of cloud droplets into raindrops 
is proportional to the average liquid water content. 

The rate of precipitation is parameterized by 

r = K(p)! f (40) 

where K(p) is an empirical function which varies with 
height. Because our model of clouds does not explicitly 
contain the vertical velocity in the clouds, it is not __ c (, 
possible to include detailed cloud physical processes.-'^ 

With these assumptions, we obtain a close 'set of .. 

•equations for the ten unknowns; M c , S, 6; h c , s c , Qc, /»_• 

■ C, r, e. The equations are as follows: 


d 

Qi-Qn = L{c— c)H [_M c (j e — 3)3 (41) 

dp 

d 

Q 2 =L(c-e)-L—tM c (q~gn (42) 

dp 



dp 

r=ld 

e=i 


(43) 

(44) 

(45) 

(46) 

(47) 

(48) 


and the saturation relations, similar to (33) and (34), 

which relate s c , q c and k c . 

In the above, only the terms Q->, s and 

q are directly observable by means of synoptic-scale 
soundings. ; 

The simplified model of the cloud ensemble is 
schematically illustrated in Fig. 3. The ensemble ex- 



Fig. 3. The idealized model of a cumulus cloud ensemble 
(see text for notation). 
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changes Its mass., heat, moisture and liquid water with 
its environment. Both entrainment and dctrainmeirL^ 

can take place at the samejieight ~ "V 

Using (43), (44) /'(T5) > ve further modify 

(41) and (42) to obtain * — 

ds 

Qi-Qs^-Me Le, (49) 

dp 


0 2 - LM C M(q* ~U) - Le. (50) 

d P 


These equations were derived by Arakawa (1971), 
Ooyama (1971) and Yanai (1971b), 

Eq. (49) shows that the apparent heat source, in 
addition to Qr, consists of an adiabatic warming due to 
a component of downward motion which compensates M c , 
and a cooling due to the re-evaporation of cloud droplets. 
Eq. (50) shows that the apparent moisture sink is due 
to a drying effect by the compensating downward motion 
and that the delrainment of i water vapor and the re- 
evaporation of detrained liquid water act as moisture 
sources for the environment. These effects were also 
discussed by Gray (1971, 1972). 

Eliminating e between (49) and (50), we obtain 

dh 

Qi-Qi-Qn - o(7i* -l)-M c —. (51) 

dp 


Elimination of c between (44) and (45) gives 

dh c - ~ dM c 

SQL* —h) - M 0 — (A e -/0— = °> 

dp dp 


(52) 


.or 


(&{7i*—Ji)—M e — -+- ( h c —h ). 

dp dp 


(52a) 


Because of the presence of the second term on the right 
of (52a), entrainment does _not nec essarily^ give. 

dydp> 0, as in the case of entrainment models of a 
. s i n n-le cloud (e.g., Levin, 1959; Simpson and Wiggcrt, 
1969; Weinstein 1970; Warner, 1970; Simpson, 1971; 
; and many others). 

Combining (51) with (52), we confirm that 

0 - - dM c 

<?, - Qi-Qu = M ( —{h c ~7i)+ (A c - A)—, 


d 

«— [M C (7/ C ~A)]. (53) 

sp 

Vertical integration of (53) gives, using (11), 

g F=M c (k c -h). (54) 

We see that (54) is identical with (22), when only one 
type of clouds is considered. 



d. The conditions at the cloud base 

The foregoing set of equations does not apply in the 
sub-cloud layer. At the cloud base M c and h, are de- 
termined by the following considerations. 

In the sub-cloud layer between pu and po (sec Fig. 4), 


-Mqb Pb 

— 1 CLOUD BASE 


0 ■' ■ ? SEA SURFACE 

s 0 lEq 

Fig. 4. Energy transfer processes between cloud 
base and sea surface. 

the heat and moisture budgets are ' 

^ ((?1 _ W =5o-— c p (T c -T)b, (55) 

l * 

— ^-Q 2 =-'M C ]s(qe—g)n~-TEo i (56) 

g g 

where subscript B refers to the cloud base, and 

Ap = p 0 -pB. Subtracting (56) from (55), we have 

&•-*>» < 57 > 

g g 

Using observed values of Q h Q-i, Qr in the sub-cloud 
layer, and Fo, (57) gives 

M C B{hc~h)a = gTo-^P (Qi - Qi -Qn)> ^ 

= gFn = fi- 

Similarity, (55) gives 

M c bc p (T c — T)ii = gSo— ^P(Qi—Qr)‘ ( 59 ) 

-/* 

If the Bowen ratio B is known, Sb = BFq/(\+B) so that 
fi and / 2 are known. Because 

= M c ii {h c — h *+A* — h) b, 

— (l+y)/2+Afcc(^* — h)u, 



where 
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we onlain 


tmd 




/i-tt+y)/* 

h*-h 


JicB~-k — fljM cB- 


(60) 

(61) 


e. Numerical method oj solution 

After specifying the cloud base conditions, the follow- 
ing set of equations is solved for each cloud layer by 
iteration: 


(tr) ^5 

~M C — =Qi-QR+Le^" 

dp 

-O) , (►) 

P>. -fi-gF/M* 

Sc $ — L ^'c h J 

1+7 


Lit' -Tn 

1+7 


8 w (k*— Q r+Mc 


m 


dk 

dp 


» 




t- 


! 


dM c 

dp 


» 


3 (»>=(»)_. 1 3M C _ 

c w— [Af« g c ] g-^V-s) 

d/> 

(O-O) 

\ = (0 7 

r^+o =c v 

/n dp 

, =(v) 

rw=ia 

<+> = 5^/ 


— r(v) 


(62) 

(63) 

(64) 

(65) 

( 66 ) 

• (67) 

( 68 ) 

(69) 

(70) 

(71) 


In the above, » is the index of iteration. 

We start with (62), taking an initial guess e (0) =0. 
Among the above ten equations, only r the liquid water 
equation (69) is a differential equation. We note that, 

at. the final converged step, M c , §, 
negative. 


I must be non- 



3. Data 

To test the computational scheme outlined in Sect ion 
2, daily vertical profiles of Q h Qi and <2 jo which are 
representative over an area, say a 300-km square, are 
needed. Ideally, these must be supplemented by daily 
values of precipitation and heat and moisture transfer 
from the ocean surface. To the authors’ knowledge, 
such an ideal data set does not exist. 

We selected the set of upper air, surface and ship 
observations of Operation Redwing, which was col- 
lected in the Marshall Islands area during the period 
from 15 April to 22 July, 1956 (joint Task Force Seven, 
1956). This special data set was used in the studies 
of Rosenthal (1960), Nitt.a (1970, 1972a, b), M. Mura- 
kami (1971, 1972), Janota (1971) and Hastenralh 
,(1972) for various purposes. During the period of 
observation, a number of disturbances passed through 
the area. 



For practical geometrical considerations, we selected 
five stations which form a pentagonal region whose 
area is ~62X10 4 km 2 (Fig. 5). In addition, surface 



observations taken aboarcPth.e. U,S,S.-Esles , .stationed^ 
near Bikini Island, we^alsoused. At the five stations, 
all the upper air soundings were made four times a 
day (0300, 0900, 1500, 2100 GMT), while all the sur- 
face and the ship observations were made four or 
eight times a day, depending on the meteorological 
elements. To insure the accuracy of lateral flux com- 
putations, all the standard and the significant level 
radio soundings and the wind aloft data were combined 
and re-edited at UCLA. The processed data are la\ er- 
averaged winds, temperature and moisture at intervals 
of 50 mb. These were used to obtain the estimates of 
and 02 . 


4 . Large-scale mass, heat and moisture budgets 

a. Procedure 

The method for obtaining the large-scale vertical 
motion a>, the apparent heat source 0i, and the ap- 
parent moisture sink 0«> is basically similar to those 
used by Yanai (1961a, 1963), Nitta (1970, 1972a) 
and Reed and Recker (1971), but a somewhat refined 
computational scheme for the horizontal integrations 
was used in this study (see Appendix). 

The area-averaged horizontal divergence and rela- 
tive vorlicity over the pentagonal region were computed 
by 



(^<j) udy— <j> vdxj, 
udx-\- <j) vdy^j , 


(72) 

(73) 


where u and v are the eastward and the northward com- 
ponents of the winds along the periphery of the penta- 
gon whose area is denoted by A . 
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The average vertical ^velocity is obtained by 



The original estimates of V ■ V were slightly corrected t 
rnake oj vanish at 100 mb. A linear correction of V 
increasing with height, was used. 





I he rate of heating due to the apparent heat source 
Qi was computed by the first law of tile thermo- 
dynamics, 


Qi 

Cp 


dT 


dt 


-fv-vr-K 



w , « as 

+(v-rv-rv-v>+ , 

Cp Op 


at 


(7$) 


The rate of equivalent heating due to the apparent 
moisture sink @ 2 was obtained by 



The corrections on V - V and cb were introduced into (75) 
and (76) so as not to violate (74). 

Because of the lack of humidity soundings above the 
250-mb level, a constant relative humidity (50%) was 
assumed between the 100- and 250-mb levels, although 
some earlier studies (Mastenbrook, 1965, 1968) sug- 
gest higher humidity values in the tropical upper 
troposphere. 

All the computations were carried out at every ob- 
servation time (390 samples) and at intervals of 50 mb. 
In this section we mention only the time averages of all 
the samples. An investigation of the variation in lime 
•of the large-scale parameters will be reported in another 
paper. 


b. Divergence , vertical velocity and vorticily 

The vertical profile of the average horizontal diver- 
gence is shown in Fig, 6. There is a deep layer of com 



luc. 6. The mean horizontal divergence. 




vergence up to the 350-mb level and a shallow layer of 
intense divergence aloft. These features are similar 
to the results obtained by Williams (1970) and Reed 
and Recker (1971) for active cloud clusters. On the 
average, there is no evidence of concentrated mass 
convergence in the boundary layer. 

The average vertical profile of the mean vertical p 
velocity is shown in lug. 7. The upward velocity 



FiO. 7. The mean vertical p velocity. 

reaches its maximum at the 350-mb level. 

As a by-product, vertical profiles of relative vorticity 
were also obtained. The time averaged profile is shown 
m Fig. S. The vorticity is positive in the lower tropo- 



Fig. 8. The mean relative vorticity. 

sphere and strongly negative in the tipper troposphere. 
*1 hesc features are similar to those computed by Yanai 
(1963) using the 1958 test. data. The profile is also very 
similar to those obtained by Williams (1970) and Reed 
and Recker (1971) for active cloud clusters. 

c. Static energy 

Mean vertical profiles of the dry static energy s, the 
moist static energy h and its saturation value h* are 
shown in Fig, 9. All the values below the 450-mb level 
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Fie. 9. The mean dry sialic energy s, moist static energy h 
(both solid), and saturation moist static energy h*' (dashed) of 
the environment. 

are systematically higher than the corresponding valu^ 
obtained by Riehl and Malleus (1958) for the mean 
equatorial trough zone, and by Hastenrath (1971) for 
Palmyra and Christmas Island. 

In the lower troposphere h is much larger than s, 
■owing to the presence of moisture; /Mias a minimum at 
-625 mb. In the upper troposphere h approaches s and 
the distinction between the two curves becomes 
negligible above the 200-mb level. The existence of such 
a minimum of It, with the upward increase of k in the 
upper troposphere, led Riehl and Malkus (1958) to 
hypothesize a process of selective upward transport of 
heat by "hot towers.” 

Note also that values of h in the lower and middle 
troposphere are much below the saturation values h* 
which show a marked constancy in the layer between 
300 and 600 mb. In contrast to the other large-scale 
parameters, the day-to-day variations of 5 and k are 
small regardless of the passage of disturbances. This 
fact suggests that cumulus convection occurs in re- 
sponse to the large-scale forcing in such a way that the 
observed profile of h is maintained. 


d. Apparent heat source and apparent moisture sink 

The time- averaged vertical profiles of (?i and Q> 
are shown in Fig. 10. Roth are expressed in units of 



F:c. 10. The mean apparent heat source Q : (solid) and moisture 
sink Q i (dashed). On the left is the. radiational heating rate given 
by Dopplick (1070). 




heating rate (°C day -1 ). Note that there is a consider- 
able difference between them. This was alreay noted 
in a study of typhoon formation (Vanai, 1961a) and 
was used by Wada (1969) and Nitta (1972a) to measure 
the activity of cumulus convection. The averaged ()j 
has a maximum (6.4 C day ' 1 ) at 475 mb. On the other 
hand, the average Q-> has its principal maximum at a 
much lower level (775 mb) and a secondary peak at 
525 mb. The profile of Qi is very similar to those ob- 
tained by Reed and Recker (1971) arid Nitta (1972a). 

With the daily vertical profiles of Q\ and Q~>, we can 
obtain the required vertical tlux of eddy heat trans- 
port, F — — (kW)/g, provided that the daily vertical 
profiles of the radiative heating rate Qn are known. 

Unfortunately, wc have no direct measurements of 
Qjt in conjunction with the upper air soundings. In 
this experimental study, we use the climatological 
vertical profile of Qn given by Dopplick (1970). The 
profile averaged from his Figs. 26 and 27 (between the 
equator and ION) is shown on the left side of Fig. 10. 
The net radial tonal cooling in the troposphere is on 
the order of 1-2C day -1 . This profile is very similar to 
the profiles of infrared cooling alone, as measured in 
the tropics with the Suomi-Kuhn infrared radiometer 
(Riehl, 1962; Cox, 1969; Cox and Hastcnrath, 1970). 

The required vertical eddy heat flux F is obtained by 
the integral (11). The flux is assumed to be zero at the 
100-mb level near the tropopause. The mean vertical 
profile of F is shown in Fig. 11 in units of cal cm -2 



Fjg. 11 . The derived vertical eddy heat flux. 

day -1 . The vertical flux has a maximum at 700 mb, a 
slight minimum near 900 mb, and its largest value (389 
cal cm -3 day -1 ) at the sea surface. 


c. Calibration of the large-scale budgets 

Although the average value of F 0 falls in the range of 
previous estimates of sensible plus latent heat flux by 
Garslang (1967) and Holland and Rasmussen (19/C 
we do not achieve good agreement with the direct 
.. surface observation of precipitation nor with an esti- 
mate of the evaporation using a bulk aerodynamic 
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(77) 


( 78 ) 



> t,ie ab °ve, P, C q q, 0 and V are the density of air in 
the surface layer, the empirical bulk aerodynamic 
coefficient for latent heat flux, specific humidity, p 0 - 
re, and the absolute air velocity. The 
subscript slOynd a refer to sea surface and anemometer 
level respectively. 

i he meteorological data of Operation Redwing have 
several deficiencies for the computation of (77) and (78) 
OnJyj?ne ship near Bikini Island, the U.S.S. Ectes (see 
5), had suitable records for the bulk aerodynamic 
computations. However, since most of the islands are 
small in size, it was assumed that by using sea surface 
temperatures extrapolated from monthly mean charts 
(U.S. Department of Interior, 1968), standard surface 
observations taken at the island stations would also 
provide a rough estimate of the evaporation and the 
Bowen ratio over the pentagonal region. All the surface 
observations (eight times per day) were used. 

The Bowen ratio was computed for each of the nearly 
800 observations. The time-averaged Bowen ratio for 
the pentagonal region was 0.076. This value agrees 
favorably with the previous estimates made for the 
eastern Caribbean by Bunker (1960), Garstang.(l%7)"' 
and Holland and Rasmusson (197$y — ~~ — 

A value of 1.23X10~ 3 was adopted for the 8 m bulk 
--dy-mic coefficient C q , as suggested by Pond el al. 
(19/1). It has an uncertainty of about 20%. The 
anemometer aboard the Es(es was mounted at 14 m 
above the sea surface so that a correction was applied 
to C g . The resulting average evaporation for the 
pentagonal region was 0.24 cm day 1 2 . The Estes. data 
gave an evaporation rate of 0.33 cm day 1 . 

Ibere are a number of ppssible reasons for the dis- 
agreement between the large-scale budgets and the 
surface observations. 

1) The radiational cooling rate given bv Dopplick 
may not be adequate for this region which 'has a large 
amount of cloud cover and active cumulus convection. 

4 he actual cooling rate is likelv to be much less. 

2) The observed surface precipitation values are 
likely to be too low. In the original data tabulation, 
there is no clear distinction between small rain amounts 
and missing observations; they had to be assumed zero. 

3) An examination shi. ws that the geographical dis- 
tribution of the computed evaporation is almost 
entirely the result of the surface wind distribution. Ac- 
cra ding to the Joint Task Force Seven (1956), the ob- 
servation sites of Eniwetok, Rogcrik, Kwajalein, Wake 
and Tarawa are excellent and entirely unobstructed. 
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rhe Majuro and Kapingaraarangi sites arc excellent 
)Ut are partially screened by palm trees. 1 onape, 
Cusaie and Truk have island obstructions. NearS v 
ill the stations with partial screening or island obstruc- 
tions showed significantly less variance in the total 
rind velocity. We believe that this island eflect re- 
sulted in an underestimate of the evaporation on the 

:>rder of 0.1 cm day -1 . . . A 

4) Despite the use of a sophisticated method en - 
ployed for the line integrations (Appendix), the 
average distance between the upper air stations was 
probably too large to obtain accurate estimates of the 
mass, heat and water vapor fluxes, ibis makes the 
estimates of Qi and Q-> inaccurate. 

We therefore conclude that our estimates of 
and Q 2 mav have errors of a factor 1.5- 2.0. But in the 
following test computations of the cumulus ensemble 
properties, we do not apply any corrections on the 
values of (>,, Qi and Q, £ . Thus, there will be an internal 
consistency of tne precipitation amount derived irom 
the model with the initially determined heal and 
moisture budgets. 

5. Average properties of Cumulus ensemble s 

a. Finite-difference computational scheme 

The equations mentioned in Sections 2d and 2e 
were applied to a 18-laver model of the cumulus en- 
semble. AU the vertical differentiations were replaced 
by finite differences. The cloud base was assumed to be 
at 950 mb. Prior to the. computation, the input data ot 
Q y and Qi were subjected to a low-pass filter to remove 
diurnal variations, which smoothed out some of he 
noise while reducing the number of samples from 590 to 

386 

In addition to the cloud base conditions (with a 
fixed value of the Bowen ratio of .0/6), we need a 
boundary condition for the liquid water equation 

cloud base. 


We applied I60s. (62)- (71) to the 386 cases. The i( 
ition converged^ 366 cases. Detailed studies have not 
7 e t been made of the individual cases. In the following 
iVe show only the time-averaged vertical distributions ol 
the 366 samples. 



5. Cloud mass jinx 

Fig. 13 shows the average vertical profiles of the cloud 
mass flux M r ., the large-scale mean mass flux M — —a>> 



and the residua) mass flux in the environment,. 

The most important feature seen m this diagram is 
that M e exceeds M excepting near the tropopausc; i.e ; , 
the upward mass flux in active cumulus clouds is 
larger than the mass flux required from the large-scale 
horizontal convergence. Another significant feature is 
that Mr. is very large near the cloud base. Ihis sug- 
gests that a large number of shallow cumulus clouds 
co-exist in the region with taller cumulonimbi. Re- 
member that M c = ^m m >' (d5)< 

As a result of the large M C} the residual mass flux is 
generally downward. This sinking motion warms the 
environment by adiabatic compression. 

c. Entraimpent and detrainment 

The average vertical distributions of entrainment and 
detrainment are shown in Fig. 14. “Because of the co- 



Fig. 14. The average detrainment 5 (solid) and 
entrainment < (dashed). 


existence of a variety of cloud types, entrainmenUo 
and detrainment 5 lake place at the same height. 'Ihe 
average entrainment is largest in the lowest layer, with 
a weak secondary maximum centered at -about 4o0 mb. 
The entrainment goes to zero at the top; of the tiopo- 

( . h 5 
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sphere. On the other hand, the average detrainmen t has 
two strong maxima, i.e., in the lowest layer and near 
the 200 -mb level. The lower maximum must be the 
result of large number of shallow clouds that detrain 
immediately above the cloud base. 


d. Moist static ejiergy of clouds 

The mean vertical profile of the moist static energy of 

the cloud ensemble, li c ~ (Si [see (37)], 

shown together with the average profiles of h and h * in 
Fig. 15. There are several interesting features in the 
shape of this profile. 




0 5 


col qnrf 1 


Fig. 15. Left: The average moist static energy of the cloud 

ensembles hc (solid), compared with the moist static energy ft 
(thin solid) and saturation moist static energy ft* (dashed) oi the 
environment. Right: Excess moist static energy of the emud 
ensembles. 


Near the cloud base, h c (84.0 cal gm J ) is almost 
equal to 7j*, the saturation value of static energy of the 
environment. Thus, the cloud air has practically zero 
buoyancy. This fact is consistent with the early ob- 
servation of the tradewind moist layer made by Malkus 
(1958) and a more recent report by Simpson and 
Wiggert (1971). ~ 

Another noteworthy feature is that the average h c 
increases with height in the upper troposphere. Because 


of the definition of h c , this indicates that wc are sam- 
pling less diluted, warmer clouds in the upper tropo- 
sphere. Near the tropopause 7i c has nearly the same value 
as at the cloud base. This tends to confirm the existence 
of undilute “hot towers” hypothesized by Riehl and 
Malkus (1958). Extending a spherical vortex model of 
Levine (1959), Malkus (1960) showed that large cloud 
elements can reach the tropopause without much 


dilution. 

On the right-hand side of Fig. 15, the e.xccss static 

energy h c ~~h is shown. The excess reaches the maximum 
(3.7 cal gmr 1 ) near 65&-700 mb. 





c. Cloud temperature, moisture and liquid loafer content 
In Fig. 16 the average profiles of the excess tempera- 
ture T c —f, the excess mixing ratio q c — q, and the liquid 
water content l of the cloud ensemble are shown. 



Fig. 16. The average excess temperature, excess mixing ratio, and 
liquid water content of the cloud ensembles. 

The excess temperature of the cloud rapidly increases 
from — 0.1 C at the cloud base to 3. 1C near 500 mb. 
Corresponding to the secondary maximum of the en- 
trainment (see Fig. 14), the excess temperature de- 
creases with height from about 500 to 250 mb, and then 
increases to the maximum near 200 mb. The comparison 
of the derived excess temperature with observations is 
difficult because there is very little literature on cloud 
temperatures excepting the classical data obtained for 
shallow trade cumuli in the Caribbean (roughly be- 
tween 800 and 930 mb), which were reported by Malleus 
(1952, 1954, 1958) and Bunker (1959). According to 
Malkus (1958), the excess temperatures of active 
buoyant part of these cumuli are on the order of 2C. 
On the other hand, there are also active downdrafts in 
the clouds, which have typical temperature deficits of 
— 1C. We feel, therefore, that the derived excess tem- 
perature is probably slightly overestimated. 

The magnitude of the excess mixing ratio (the maxi- 
mum value is 5.4 gm kg -1 at 750 mb) also seems to be 
too large when compared with the previous observations 
(Malkus, 1954, 1958; Bunker, 1959). However, Bunker 
d al. (1949) reported a sounding made in a cumulus 
cloud, which showed an excess of 5 gm kg -1 at about 
1.4 km height. From a balance requirement for the 
vertical transport of water vapor in the tradewind 
cloud layer, Malkus (1962) estimated an excess mixing 
ratio of 5.6 gm kg -1 at the same height. 

The average liquid water content (gnvm'" 3 ) of the 
cloud ensemble looks very reasonable, despite our crude 
parameterization of the precipitation process. There are 
a rather large number of studies on the liquid water 
content of maritime cumuli (Warner, 1955, 1970; 
Squires, 1958; Warner and Squires, 1958; Ackerman, 
1959; Telford and Warner, 1962; Simpson el al., 1965; 
Simpson and Wiggart, 1969). Ackerman (1963) mea- 
sured the liquid water content in hurricanes. 

It is customary to discuss the vertical distribution of 
the relative liquid water content (the ratio of mean liquid 
water content sampled in clouds to its adiabatic value) 
as shown in Fig. 17. Wc note that the computed rela- 
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Fig. 17. The average relative liquid water content (l/l a ) of the 
cloud ensembles (computed value, thick line), observed results 
(thin lines). 

live water contents in lower levels are much higher than 
the observed values in shallow cumuli reported by 
Warner (1955), Squires (1958) and Ackerman (1959). 
The highest altitude of observation made in hurricanes 
by Ackerman (1968) is 18,000 ft above the cloud base. 
It is interesting that the high-level portion of our com- 
puted curve seems to join with Ackerman’s hurricane 
curve rather smoothly. 


/. Condensation, precipitation and re-evaporation 

It was shown in Ecjs. (46) and (47) that the liquid 
water equation for the cloud ensemble is 

_ dM j - _ 

~o7+: +C=K(p)l f 

dp - 



in which the rate of precipitation (the right-hand side) 
is very crudely parameterized. Formally, the form 
resembles those used by Kessler (1965, 1969) and 
Ogura and Takahashi (1971). However, in our model 
this term gives the rate of precipitation per unit mass 
of air (resulting from all clouds) as a function of 
pressure, so that a number of factors such as the cloud 
7 ^uiss41ia-x should be implicitly involved in K(f>). After 
K • several experiments to Lest the convergence of the 
scheme, K(p) was assumed to increase nearly expon- 
entially from 0 at. the cloud base to 60.0 day -1 (~10~ 3 
sec -1 ) in the top layer. The assumption is admittedly' 
crude and empirical. However, this crude treatment 
still yields rather reasonable results for the vertical 
distribution of the generation of raindrops. 

In Fig. 18, the vertical distributions of the mean rate 
of condensation, precipitation and evaporation of clouds 



Fig. 18. The average rates of condensation Qr- (solid), precipita- 
tion Q„ (dashed) and re-evaporation Qk (dash-dotted), in equiv- 
alent heating units. 


are shown together in equivalent heating units. As 
required by the large-scale heat budget, the release 
heat associated with the precipitation is nearly 1000 
cal cm -2 day r_1 . In the upper troposphere the curve of 
precipitation is shifted upward by the vertical transport 
of liquid water. In the lower troposphere the rate of 
generation of raindrops is much smaller than the rate 
of condensation. Instead, a large amount: of re-evapora- 
tion of clouds is taking place. The efficiency of the rain 
generation process (the ratio of the total precipitation 
to the total condensation) is about 84%. 


3 //V 


r 


g. Heal balance of ike environment 


Having obtained all the variables which characterize 
the cloud cluster, we can discuss the observed apparent 
heat source in terms of the interaction between the 
clouds and the environment. As shown in (49) the 
heat source is expressed as 


ds 

Qi-Qr=-M c Le, 

dp 

be., the sum of the adiabatic heating due to the com- 
pensating downward motion in the environment and 
the cooling due to evaporation of cloud droplets. These 
are also discussed by Gray (1971, 1972). By assumption 
(38), there is no detrainmenl of sensible heat from the 
clouds to the environment, except in the uppermost 
layer. These three features are shown in Fig. 19. 
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Fig. 19. The observed heat source Qi—Qir (dash-dolled), 

. adiabatic heating by the compensating downward mass flux 
— Af,;(ds /dp) (solid), evaporative cooling —Qk — ~*-e (dasnecj, 
and detrain men l of heat 5(7 0 — ■ T) (thin solid). 

The clouds act as a beating agent primarily through 
Ahcir induction of Ike compensating descent ( — M c ). This 
interpretation docs not mean that the actual mass flux 
in the environment, il I~M — M c , must necessarily be 
downward. However, nearly all the computations show 
that M is negative. 

The second important contribution to the apparent 
heat source is the cooling due to re-evaporation of the 
detrained cloud droplets. This effecl is pronounced in 
.the lower troposphere. 


i 




:h. Moisture balance of ike environment 

The average moisture balance of the environment 
Is shown in big. 20, in units of gm kg -1 da}- l , and 



dg . d °l 

+V- (gV)+— , 
dt dp 


dq - 

= -K c — M(2*-y+f)< 

dp 

The apparent moisture sink (Q-i) is primarily due to 
the induced downward motion ( M c ). This factor 
alone would drv out the environmental air. It is the 
detrainmen t of water vapor and of liquid water from 
the clouds which makes up the balance needed to 
maintain the moisture of the environmental air. 

This would imply that the shallower, non-precipitat- 
ing cumulus clouds are needed to help the growth of 
the deep, precipitating cumulus towers by supplying 
the mejsiture. It is known from photographic observa- 
'? tions that large-sized clouds are surrounded by great 
numbers of small cloudlets (Malkus, 1958). 



6. Conclusions and future work 

We have shown that by combining the large-scale 
heal and moisture budgets and an entramment- 
detrainment model of a cumulus c ourt eirsemb e < 
large amount of information can be obtained about the 
bulk properties of tropical cloud clusters. 

Because of tbc uncertainty inherent m the data used 
for our illustrative example, the derived quantities, 
such as the cloud mass flux, the excess cloud temper - 
lure, etc.., in this particular case, may be m erroi by a 

^Tbcmost. important conclusion of this study is that 
the cloud moss flux exceeds the mean vertical mass flux 
required by the large-scale convergence, Urns causing a 
. compensating sinking motion between active Ms. Me 
large-scale healing of the environmental air is 
due to its adiabatic compression m the com pm sal; g 
downward motion. The cooling due to re-evaporal ion f 
liquid water detrained from the clouds is also an importe d 
factor inike heal balance of the environment. Counteracting 
the drying due to the environmental sinking motion are the 
large amounts of water vapor and liquid water which a 
detrained from the clouds, especially from the shallow 
clouds in the lower troposphere. Thus, t « fallow 
non-precipitating cumulus clouds, by the 
transfer of water vapor, support the grow h o the . de 1 , 
precipitating cumulus lowers. Some of the above con 
elusions were also reached by Gray (19/2) from an 

"uX '““discussed only the time-averaged 
vertical distributions of the cloud cluster parameter. 

An examination of the day-to-day changes of these 
parameters together ■ whh the » 

between cumulus 

convection and the large-scale tropical motions- , 

In this section, we wish to point out several de forts 
in the reported study, which arc subject to future 
improvements. 

CO The chosen pentagonal region may be too large 
to he viewed as a portion of large-scale motion systems. 
As evidenced by the analyses of rainfal and evapora- 
tion amounts,' the region contains substantial hori- 
zontal variations of the large-scale pmamjgers In 
particular, the possible horizontal gradient of the lax£ 
scalc vertical motion causes an ambiguity in tl 
pretation of /iV. We plan to sub-divide the region i to 
reveal the regional differences ol the cumulus cluster 
properties between the 1TCZ and the relatively dry 
tradewind regime. 

(ii) There is some evidence which suggests 111 
islence of short-lived mesoscale organization of cumu- 
lus convection (e.g„ Zipser, 1969) If this is true m 
general, further modification of the formulation cull be 

ne (?h) d The weakest part of the present study is the 
crude parameterization of precipitation processes Un- 
less the area ratio a is specified we have no wa> o 
tain the vertical velocity in the cumulus clou « 
hence the detailed cloud physical processes which were 
bypassed in this study. 
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Finally, we stress the need for better, carefully 
designed observational programs. To apply our method 
of computation, better data sets are needed to obtain 
accurate estimates of the heat and moisture budgets 
over, say, a 300-km square area. This must be accom- 
panied by equally accurate measurements of radia- 
tiona! cooling. The measurements of precipitation and 
evaporation also needed to be improved, and the in- 
ternal structure of the cloud dusters should be directly 
sampled by aircraft. We hope that the currently 
planned GATE experiment will yield such a data set. 
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APPENDIX 

Objective Analysis of Large-Scale Field Variables 

Stations Eniwetok, Ponape, Kusaie, Majuro and 
Rongerik form a pentagon (Fig. 5). Let x if yu (i= 1, 2, 
5) be the longitudinal and latitudinal distances 
relative to the center of gravity of the pentagon. 

Let Hi (?— 1, 5) be an arbitrary scalar field 

variable (zonal wind component, for example) observed 
at the stations. 

We assume that 

Ui = A Xi 2 + BxijiA-Cy r+ Ey { -f F i 

(i= 1,2, ---, 5). (Al) 

Because we have six coefficients to be determined 
from five observed values of we impose the additional 
constraint that the curvature of the quadratic surface 
should be minimized as in the case of skilled hand 
analysis. Thus, we minimize 


/ d 2 u \ 2 fdht dhi\ 2 

(W + fc-i (A2) 

Eq. (Al) may then be written 

A xf+Bxiyrt-Cyf+Dxi+Eyi = u t - F. (A3) 

The imposed condition is expressed as 
SB /3A 0C\ 


dt> /a A d C\ 

^ J=o. 

OF \dF dF/ 


(A4) 



Solving 


c U OB dC dD 

— xFA — v,;yH yd A Xi+ 

OF dF dF OF 


OE 



(AS) 

we obtain 

dA/dF , dB/dF and dC/dF . 

Because 




04-C)/3C cM 


\aF aF, 


— — =^4-C), 


(A6) 


OF 

(Al) reduces to 

A (xp+bxiyi) A~C (y, 2 — bx/yi) -f Dx, FEyrbF 

(^ 1 , 2 ,..: 



5). (A7) 


Since (A7) has only live unknown coefficients, it is 
solvable from five observed values of 

All the field variables reported .in Section 4 were 
anal) r zed by this method. 

Line integrals such as f udx, $ u<iy, and the areal 
average «— (‘xrezr v )J'J'u<1xdy were exactly obtained 
by using the analytical expression (A7). 

The results of the mass flux computation (yfiudy 
— jF vdx) using the present method were compared with 
those obtained with the conventional linear interpola- 
tion method. For the time -averaged flux, both methods 
give almost identical results, but the present, method 
gives smoother results for the daily values. 
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